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1 Introduction 

The correlation between CR and LCC variations, which led to the introduction of 
a new scientific subject - 'cosmoclimatology', was found more than 10 years ago 
( Svensmark and Friis-Christensen, 1997; Palle Bago and Butler, 2000; Svensmark, 
2007 ) . The proponents of the causal connection between CR and LCC point out 
a number of facts. Firstly, there is the positive character of the correlation, i. e. 
an increase of the CR intensity is accompanied by an increase of LCC and vice 
versa. Secondly, the peak to peak amplitude of the Global LCC variations (~2%) is 
much higher than the amplitude of the variation of the energy flux, delivered by the 
Sun, or the sun's irradiance (SI,~0.1% peak to peak over the Solar Cycle), which 
requires us to find the mechanism which would explain such a large magnification. 
Thirdly, there is an effect in the LCC similar to the latitude effect in CR, i.e. LCC 
variations during the 11-year cycle of solar activity are less in the tropics than at 
higher latitudes. The same reduction of variations in the equatorial regions exists 
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also in CR due to the higher geomagnetic rigidity near the equator. 

The opponents of the causal connection between CR and LCC have put forward 
different arguments. Firstly, the positive correlation of CR and Cloud Cover (CC) 
is noticed only for the LCC, i.e. for clouds below 3 km above sea level. No signif- 
icant positive correlation has been found for higher clouds. Secondly, there is an 
altitude dependence of CC and CR, but it changes sign. If one thinks about the 
ionization of the air as the mechanism of the CR influence on CC formation ( the 
usual assumption ), then the maximum of the CR flux and ionization is at heights 
of ~12-15 km and not below 3 km, where the effect is claimed. Thirdly, there were 
no changes of CC noticed after the significant release of radioactivity during the 
Chernobyl disaster or during ground-based tests of nuclear weapons ( Erlykin et al., 
2009 ). Fourthly, there were no CC changes found during and after strong short-term 
variations of the CR intensity ( Forbush decreases or GLE - ground level events ) 
( Kristjansson et al., 2008, Sloan and Wolfendale, 2008 ). 

The purpose of the present paper is a further analysis of the possible origin of LCC 
and CR correlations found in the work of Svensmark and Friis-Christensen, (1997) 
and Palle Bago and Butler, (2000). 

2 Input data 

As input data on the CC we used the same observations by meteo-satellites incor- 
porated in the ISCCP program ( ISCCP, 1996 ), which were used in Svensmark and 
Friis-Christensen, (1997); Palle Bago and Butler, (2000); Svensmark, (2007). We 
analysed monthly means for the fraction of the total observed area occupied by the 
clouds (D2). Following the classification of the cloud heights adopted in the ISCCP 
they were classified according to the pressure at their top border as: low ( LCC, 
>680 hPa ), medium ( MCC, 440-680 hPa ) and high ( HCC, <440 hPa ). Due to 
the continuing dispute on the quality of ISCCP radiometer calibrations after 1996 
( Marsh and Svensmark, 2003 ) we started the analysis using data obtained only 
during the 22nd cycle of the solar activity ( July 1986 - December 1995 ), but later 
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added also the 23rd cycle and used the whole set of data available from then on. For 
the comparison with CR variations we used as a proxy of the Global CR intensity 
just the neutron counting rate of the Climax neutron monitor, situated at a latitude 
of 39. 4° iV ( WDC neutron data ). The CR variations at other latitudes, though 
having different amplitudes, have the same temporal behavior. The differences in 
amplitudes of the variations do not influence the value of the correlation coefficient. 

In the analysis of the latitude dependence of CR and CC variations the entire latitude 
range from -90° to 90° was divided into 9 equal intervals of 20° width. We analysed 
also the temporal behavior of the Global CC, i.e. averaged over the Globe. For 
the more distinct revelation of the non-trivial variations of CC in most cases we 
subtracted seasonal variations of CC from winter to summer, but in special cases we 
analysed also total variations including seasonal ones. Seasonal CC variations were 
calculated as deviations of the monthly mean CC values in the D2 series from the 
yearly mean values averaged over all similar months ( January through December ) 
used in the analysis. 

3 Results 

3. 1 The altitude dependence of the Cloud Cover 

The mean values of the Global CC during the 22nd solar cycle are (28.09±1.06)% 
for LCC, (19.52±1.70)% for MCC and (13.35±0.61)% for HCC. One can see that 
CC goes down with increasing altitude, which is opposite to the rising behavior of 
CR, see e.g. Hayakawa, (1965). Numerous models have been proposed to explain 
this different altitude dependence and justify the causal CR-CC connection ( see eg. 
the bibliography in Kirkby, 2007 or the recent paper by Kudryavtsev and Yungner, 
2009 ), but in our view there is still no convincing proof of their validity. The different 
altitude dependence of CC and CR is a problem for the concept of causal connection 
between them and requires further study. 

The most likely part of CR which can be connected with cloud formation is their 
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charged component, which produces ionization and which could in principle give rise 
to the growth of condensation nuclei. Balloon studies of temporal variations of the 
charged CR component at different atmospheric altitudes show that the correlation 
between variations of the charged particle flux and the counting rate of ground-based 
neutron monitors, which is rather high in the stratosphere above 15 km, decreases 
below 6 km ( Bazilevskaya et al., 2007; Bazilevskaya et al., 2008; Ermakov et al., 
1997 ). The correlation coefficient at altitudes below 3 km becomes as low as ~ 0.2. 
So it is hard to expect that CR variations observed with neutron monitors could be 
the cause of LCC variations via ion production. 

3.2 The time lag between CR and LCC temporal variations 

In Figure 1 the temporal behavior of the CR intensity (a) and Global LCC (b) are 
shown for solar cycle 22 ( 1986 - 1996 ). 

Figure 1. 

As an illustration of the CR behavior we have taken the data of the Climax neutron 
monitor. The qualitative correlation between CR and LCC can be seen by the naked 
eye: both of them reach their minimum at about the same time - July-October 1990. 
One can ask whether it is possible to find the best fit time lag between these curves, 
for which the least-squares y 2 between them has a minimum. One can imagine that 
if, say, CR variations start after the LCC ones, then CR can hardly be the cause of 
the LCC variations. The value of x 2 was calculated as 

^ Af )- E -O ii + Ai) -(lf) (f,))2 (1) 

Here LCC(U + At) is the Global low cloud cover value at the time U + At , (LCC) 
is its mean value in the studied time interval 1986-1996, Iciiiti) and (Icr) are the 
CR intensity and its mean value respectively, ndf (number of degrees of freedom) is 
the number of months taken in the analysis, At is the time lag between LCC and 
Icr, for which we find the minimum of x 2 . 
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In Figure lc we show the value of x 2 / n df as a function of this time lag within 
a At = ±1 year time interval. It is seen that x*/ndf has a very flat and broad 
minimum within -11/ +6 months time lag and it is not possible to say which of CR 
and LCC variations start first. 

3.3 Long-term variations and the fraction of LCC which correlates with CR 

As long-term variations we call deviations of LCC and CR values from their means 
obtained by averaging over the entire analysed time interval ( the dotted lines in 
Figures la,b ). In Figure 2a we show the correlation plot for the variations of LCC 
and CR during the 22nd solar cycle. 

Figure 2 

It is seen that deviations from the mean values of LCC and CR correlate positively 
with each other. The slope of the linear regression line is 0.157±0.023 and the cor- 
relation coefficient is 0.538±0.047, which confirms the positive correlation between 
LCC and CR, found in Svensmark and Friis-Christensen (1997), Palle Bago and 
Butler (2000). 

If it is assumed that CR are responsible for just a fraction of the LCC and they 
are the only agent creating this fraction, then from the observed correlation it is 
possible to estimate how large this fraction is. Such an estimate depends on the 
model of the connection between CR and LCC. Let us assume that this connection 
can be fitted as 

Y = a + bX c (2) 

where Y = ^cc) an< ^ ^ = (7ck> ' ^" ere LCC and Icr are values of LCC and 
CR intensity and (LCC) and (Icr) are their mean values respectively. The first 
and second terms in this expression determine the parts of LCC independent and 
dependent on CR respectively. If the connection between CR and LCC is linear, 
i.e. c = 1, then the slope of the linear regression line, b = 0.157, gives the fraction 
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of LCC connected with CR as ~16% as the best estimate and which should not 
exceed 20% at the level of 2 standard deviations. However, the determination of 
the a, b, c coefficients by the least-squares method shows that the best-fit connection 
between CR and LCC is non-linear rather than linear. The derived values are a = 
0.9783 ± 0.0008, b = 0.0177 ± 0.0010 and c = 8.65 ± 0.46 ( full line in Figure 2a ). 
This shows that the most likely fraction of LCC connected with CR, which can be 
derived from expression (3), does not exceed 2% around X — 1. 

This conclusion is valid only if the models of the CR and LCC connection are true 
and CR variations at the Climax latitude of ~40°N are a good representation of the 
Global CR variations. For values of c < 1 the fraction of LCC, which varies together 
with CR, can be higher. Unfortunately, due to the relatively small magnitude of 
the CR and LCC variations, it is impossible to distinguish between the models of 
the connection from the LCC-CR correlation plot of Figure 2a. Although the least- 
squares method gives preference to the value of c > 1, in the region where there are 
experimental data the behavior of curves for different values of c and corresponding 
least-squares sums differ insignificantly from each other. 

The two most popular models adopted for the connection between LCC and CR, 
which are discussed in the literature, will be considered. They are based on the 
connection between the ionization rate q and ion density n in the atmosphere. The 
first model assumes that n oc y/q, the second one - n oc q ( Mason, 1971; Bazilevskaya 
et al., 2008 ). Applied to the LCC-CR connection they correspond to c = 0.5 for the 
first model and c = 1 for the second one. It is appreciated that elsewhere ( Sloan 
and Wolfendale, 2008 ) the n oc y/q model was used. Were that to be adopted here 
the upper limit to the CR fraction would go up by ~ factor 2, to 40%. Conversely, 
if n oc q, the Sloan and Wolfendale limit ( Sloan and Wolfendale, 2008 ) would fall 
to 12% at the 95% confidence level. 

Experimental data for the charged CR and ion density, which is, of course, relevant 
here, give preference to c = 1 and show no evidence for a change with altitude at 
least for altitudes about 7-30 km above sea level. At lower altitudes they indicate 
the trend to c > 1, which qualitatively agrees with our best fit value of c = 8.65 
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( Ermakov et al., 1997 ). Keeping in mind all the necessary reservations we persist 
with our estimate of the fraction / < 20% for the latter model since it is based on 
the experimental data. 

The authors referred to above ( Bazilevskaya et al., 2007 ) also stressed that varia- 
tions observed with ground-based neutron monitors correlate well with charged CR 
fluxes only at altitudes above 15 km, thus, they may be correctly used as a proxy of 
ionizing component only for stratospheric altitudes. At altitudes below 3 km the cor- 
relation coefficient falls to about 0.2. Therefore it is unlikely that variations observed 
with neutron monitors can be followed by similar variations of the ionizing compo- 
nent at low altitudes. The good positive correlation between the counting rate in 
neutron monitors and LCC found in Svensmark and Friis-Christensen (1997); Palle 
Bago and Butler (2000) should have a cause different from the ionization of the air 
by CR with subsequent formation of cloud droplets on these ions. 

3.4 Short-term variations 

Figure 2a and the analysis made in the previous subsection are relevant to the 
total CR and LCC variations about their mean value, the main contribution to 
which is given by the long-term variations, connected with the 11-year cycle of 
solar activity. In order to reveal the possible correlation of short-term CR and LCC 
variations we removed the contribution of long-term variations. For that purpose 
the temporal behavior of CR and LCC were approximated by a 5-degree polynomial 
fit ( dashed lines in Figures la and lb ) and deviations from this fit were calculated. 
Since we used the D2-set, i.e. monthly averaged data, this analysis relates to the 
variations of monthly duration. We did not find any significant correlation between 
CR and Global LCC ( Figure 2b ). The slope of the linear regression line was 
b = -0.060 ±0.062 and the correlation coefficient r = -0.104 ±0.092. This negative 
result is in indirect agreement with the absence of even shorter daily-long variations 
of LCC during Forbush decreases or GLE, analysed in Kristjansson et al. (2008); 
Sloan and Wolfendale (2008). 
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The preliminary conclusion which can be drawn from the analysis so far is the 
following: if CR are resposible for a part of the LCC then it is most likely that 
this part is small, viz. less than about 20%. The absence of short-term correlations 
between CR and LCC indicates that the assumed causal connection between them 
could be revealed only on a longer time scale, not less than several months, which 
could be understood if the Global LCC has a monthly or longer inertia. 

3.5 The anticorrelation between LCC and CC at higher altitudes: MCC and HCC 

A significant argument against the causal connection between CR and LCC is the 
anticorrelation of LCC and CC at higher altitudes: MCC and HCC. In Svensmark 
and Friis-Christensen (1997); Palle Bago and Butler (2000) the authors claim that 
they cannot find any positive correlation between CR, MCC and HCC, similar to 
that found for LCC. It is true, since both MCC and HCC anticorrelate with LCC. 
This is illustrated in Figure 3 for the same 22nd solar cycle. 

Figure 3 

The left set of panels shows the temporal behavior of the Global MCC and LCC 
together with the correlation between their variations. The anticorrelation is clearly 
seen both in long-term and short-term variations. The existence of short-term anti- 
correlations between MCC and LCC is a remarkable difference with the case of CR 
and LCC. It points to a strong connection between clouds at adjacent altitudes. 

The right hand set of panels shows the same characteristics for clouds at the adjacent 
altitudes: HCC and MCC. The positive long-term correlation between them proves 
the existence of an anticorrelation between HCC and LCC. Short-term correlations 
between HCC and MCC are absent. A relevant point concerns the role of updrafts 
which play such a key role in cloud formation; they can in principle cause the LCC 
and MCC anticorrelations ( see later ). 

Turning to HCC with mean temperatures below ~ — 30°C at all latitudes, ice crystals 
are important and the Physics is different. The lack of an HCC-MCC correlation is 
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not surprising. 

The anticorrelation between LCC and CC at higher altitudes gives a strong argument 
against the causal connection between CR and LCC. It is difficult to imagine that, 
say, the rise of CR intensity could raise LCC below 3 km, but reduce MCC above 
this altitude and vice versa. 

3. 6 Seasonal variations 

In all previous figures well understood seasonal variations of CC have been removed. 
However, they can also be used to clarify the interaction of clouds at different atmo- 
spheric altitudes. Figure 4 shows the temporal behavior of MCC and LCC during 
the last decade of the century. It is clearly seen that the minima of MCC corre- 
spond to maxima of LCC and vice versa. Therefore, the anticorrelation of long-term 
decadal and short-term monthly variations of MCC and LCC, illustrated in previous 
subsections, is strongly confirmed on the intermediate yearly time scale. 

Figure 4. 

It has been argued that the strong seasonal periodicity of CC is most likely caused by 
the seasonal variation of the surface temperature T. Despite the fact that the seasons 
are opposite in the northern and southern hemispheres the surface temperature 
averaged over the Globe still depends on the season. The amplitude of the variation 
is about 3.8°C ( Figure 5 ) and its maximum is in July ( Global Surface Temperature 
Anomalies ). 

Figure 5. 

Figure 4 shows that the maximum of Global LCC is also in the middle of the year. 
If LCC and MCC variations are directly connected with variations of the Global 
surface temperature it means that either the LCC-T correlation is positive and 
there is no time lag between them longer than 2-3 months. If the thermal inertia of 
the Earth's surface causes longer time lags of about 6 months or longer the LCC-T 
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correlation should be negative. The latter possibility seems to us more likely since 
it is confirmed by the existence of the time lags between T and LCC long-term 
variations. ( see later ). The possible mechanism of the T and CC connection should 
also give an opposite sign for LCC and MCC variations. 

3. 7 CR and CC correlations in the 22nd and 23rd solar cycles 

At the beginning of the 21st century it was noticed that the positive correlation 
between CR and LCC ( Marsh and Svensmark, 2003; Usoskin et al., 2004 ) decreased. 
While the minimum CR intensity at the 23rd solar cycle around 2001-2003 increased 
compared with the previous minimum in 1990-1992, the LCC in the 23rd solar cycle 
was definitely lower than in the 22nd cycle ( Figure 6 ). 

Figure 6. 

It should be added that many studies have found cycle 23 to be 'anomalous' in a 
number of ways: lower sunspot number, but nearly no change in total solar irra- 
diance, extended solar minimum, a hump in the neutron monitor counting rate in 
2004 and 2005 etc. 

When we include the 23rd cycle into our analysis the correlation coefficient falls 
from r = 0.538 down to r = 0.390. Some authors ( Marsh and Svensmark, 2003; 
Usoskin et al., 2004 ) explained this fact by a fault in the calibration of the ISCCP 
radiometers, which occured at about 1995, and caused the continuous decreasing 
trend in the derived LCC values. We think, however, that this trend is not an 
artifact connected with that fault, but is a real physical effect connected with the 
rising surface temperature. In Figure 7 the temporal behavior of the Global LCC, 
MCC and HCC are shown together with the rising temperature during the 22nd 
and 23rd cycles of solar activity. 

Figure 7. 

With the addition of the 23rd solar cycle, the anticorrelation of LCC with MCC and 
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HCC becomes even stronger. The long-term anticorrelation between LCC and MCC 
increases from r = —0.636 up to r = —0.873, the short-term anticorrelation - from 
r = —0.510 up to r = —0.585. We argue that the strong short-term and seasonal 
anticorrelations between LCC and MCC are a serious argument against an artificial 
origin of the long-term decrease of LCC. The assumed calibration fault of ISCCP 
radiometers can hardly have a monthly occurence and seasonal periodicity. In what 
follows we shall analyse the CC behavior mostly taking into account the 22nd and 
23rd solar cycles together. 

3. 8 The latitude dependence of CC properties 

Both the CR intensity and the surface temperature depend on the latitude. In this 
connection it is also reasonable to analyse the variation of the CC characteristics 
with latitude. We show some of them in Figure 8. 

Figure 8. 

Figure 8a shows the latitude dependence of LCC, MCC and HCC. It is seen that 
there is a small minimum for LCC in the equatorial region, which could in principle 
be connected with the reduction of the CR intensity, but it is not confirmed by the 
local maxima in MCC and HCC. In the Polar regions, where the CR intensity is 
highest, there is an opposite decrease of LCC, which apparently is connected with 
the dominant influence of the atmospheric conditions, eg. low temperatures. The 
highest LCC is in the southern latitude bands with the largest part of the area 
occupied by oceans, i.e. with a relatively large density of water vapor. 

The altitude dependence of CC also does not correspond to the altitude dependence 
of the CR intensity. In most latitude bands MCC and HCC are smaller than LCC, 
which is opposite to CR with their intensity rising with altitude. All this shows that 
even if there is a causal connection between CR and LCC, its character is more 
complicated than the direct and positive connection. 

We have already mentioned in §3.3 that the Global LCC - CR correlation is positive: 
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r = 0.538. Figure 8b shows the latitude dependence of the CC - CR correlation 
coefficient. For a CR proxy we used just the neutron counting rate at Climax. In 
spite of the latitude dependence of the CR variation amplitude, the value of the LCC 
-CR correlation coefficient does not depend on the latitude due to the similarity of 
the temporal behavior of CR variations at different latitude bands. It is remarkable 
that, in most latitude bands, MCC and HCC have negative correlations with CR, in 
opposition to the positive LCC - CR correlation, which was the main argument for 
the claimed causal CR - CC connection ( Svensmark and Friis-Christensen, 1997; 
Palle Bago and Butler, 2000 ). Furthermore, the general similarity of the MCC-CR 
and HCC-CR correlations does not fit in with the idea that CR-induced ions cause 
cloud droplets because in HCC ice crystals dominate, where, as remarked already, 
the Physics is different. 

Figure 8c shows the latitude dependence of the sensitivity and correlation between 
MCC and LCC. The sensitivity of one variable to another, according to the definition 
( Uchaikin and Ryzhov, 1988 ), is the derivative of the first variable on the second in 
log-coordinates. In our case the sensitivity is the slope of the linear regression line in 
the MCC-LCC plot. One can notice two features: (i) the sensitivity of MCC to LCC 
and MCC - LCC correlation coefficient are negative at nearly all latitudes, which 
is another support of their Global anticorrelation. The negative sensitivity of MCC 
to LCC is difficult to explain in the framework of the causal connection between 
CC and CR, since the rise of the CR intensity should change CC similarly at all 
altitudes ; (ii) the highest negative sensitivity and the correlation between MCC and 
LCC is observed in tropical and subtropical regions: £ = —30°/ + 30° as well as in 
the southern latitude bands with the highest fraction of water: £ = —65°/ — 45°. 

4 Discussion 

The vast bibliography of the works which have been devoted to the problem of 
the possible connecion between CR, CC and climate is given in the comprehensive 
survey by Kirkby (2007). 
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The analysis made in the present work, as well as arguments presented in our previ- 
ous publication ( Sloan and Wolfendale, 2008 ), gives sufficient basis to argue that 
CR are not the dominant factor in the formation of clouds. Long-term, short-term 
and seasonal anticorrelations of LCC and MCC, which are strongest in tropical and 
subtropical regions, as well as in regions mainly occupied by oceans, allow us to re- 
turn to the traditional scenario of the main cause of the cloud variation, connected 
with variations of the surface temperature, humidity and wind velocity. However, 
one can ask whether the temperature or CC variations start first. 

We try to answer this question analysing the time lag between the surface tempera- 
ture Tand LCC, since the lowest heights of the atmosphere are closest to the Earth's 
surface and the LCC is most sensitive to the surface temperature. For this purpose 
we come back to Figure 7 which shows the temporal behavior of the Global surface 
temperature (panel a) and LCC (panel d) fitted by linear and 5-degree polynomial 
approximations. 

It is seen that besides its long-term rising trend the surface temperature has also an 
oscillating behavior similar to that of the LCC. Its amplitude is about 0.1°C and 
the phase anticorrelates with the phase of LCC with a considerable time lag. It is 
difficult to estimate the magnitude of this time lag because it depends on the degree 
of the polynomial fit. A more detailed analysis of the correlation between LCC and 
surface temperature as a function of the time lag shows that the minimum negative 
correlation coefficient is for the time lag of 5 months, but the minimum is rather 
broad. Since temperature variations are ahead of LCC variations, one can conclude 
that the former could be the cause of the latter, but not vice versa. 

The long-term oscillations of temperature of the order 0.1°C are observed in much 
longer time intervals ( Haigh, 2007; ACRIM ). They are usually associated with 
oscillations of the total solar irradiance (TSI) which has an 11-12 year periodicity. 
We have analysed the frequency spectrum of temperature variations for the 1880 - 
2008 time interval with the result shown in Figure 9. 

Figure 9. 
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One can notice the small peak at 0.007 month -1 frequency, which corresponds to 
~ 11-year period, coincident with the 11.87-year period of the solar cycle ( Sturrock, 
2008 ). The small amplitude of the peak and its corresponding low confidence level 
(2.1 standard deviations ) is presumably determined by the small amplitudes of TSI 
variations ( 1.7 Wm -2 ) and of the corresponding solar forcing ( 0.3 Win -2 ) together 
with the spread in '11 year' periods. We remark that the similar peak in the spectrum 
of land temperature variations is higher by a factor of 2, which is reasonable since 
the land is more sensitive to TSI ( see Figure 5 ). Interestingly, there is another peak 
in the frequence spectrum at 0.004 month -1 , which corresponds to a 21-year solar 
cycle and has much higher confidence level ( 7.1 standard deviations ). Therefore, 
periodic variations of the surface temperature and corresponding variations of LCC 
are most likely of solar or perhaps geomagnetic origin rather than CR, because the 
21-year variation in the CR rate is small. 

The most likely cause of the anticorrelations between LCC and MCC is the variation 
of convection flows of the air with temperature. The rise of the surface temperature 
gives rise to the growing temperature in the lower atmosphere and an upward con- 
vection flows with the corresponding rise of mean cloud heights. Since clouds in the 
ISCCP experiment are classified by the height of their upper borders the increase of 
their heights leads to the redistribution along their altitudes: some low clouds cross 
the 3 km border and become medium clouds. As a result, LCC decreases with rising 
temperature. This trend is clearly seen in Figure 7, where the rise of the Global 
temperature in the (a)-panel is accompanied by the fall of LCC in panel (d). 

The heating of the atmosphere is a slow process. The slow updraft of the air is 
the cause of the time lag between the variations of the surface temperature and 
CC. Perhaps this time lag is the reason why the maximum of the seasonal Global 
temperature, which is in July ( Figure 5 ), corresponds to the maximum in the 
seasonal variations of LCC, which is also in the middle of the year ( Figure 4 ), in 
spite of their anticorrelation in the long-term scale. The slow, but steady, upward 
convection flow of the heated air is the possible mechanism for the magnification 
of small TSI variations (~ 0.1%) giving rise to larger variations of cloud heights 
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(~ 0.7%). The slow fall of the low cloud top pressure is seen in Figure 10. The 
slope of the linear fit b is definitely negative: b = —0.148 ± 0.081 hPa-year -1 . It 
corresponds to a rise of the mean low cloud top height by about 40m in 20 years. It 
is very slow but one should keep in mind that the observed mean cloud top height 
(~2.6km at 730hPa) is only ~600m below the border (~3.2km at 680mPa) between 
LCC and MCC and any change in the height of clouds can cause the variation of 
the magnitude of LCC and MCC. 

We argue that the positive correlation of CR and LCC found in Svensmark and 
Friis-Christensen, (1997) and Palle Bago and Butler, (2000) is not evidence for a 
causal connection between them, but the consequence of a parallel influence of the 
common source - the solar activity on CR from one side and CC the other. 

Concerning the relationship between CC and ground level temperature changes, 
there have been a number of studies for particular regions. Data for the USA covering 
the period 1900 to 1990 ( Barry and Chorley, 1998 ) give roughly a total 'mean 
annual cloud cover' change of -1.5% over the 11-year cycle with an associated 0.2°C 
change in ground level temperature. A value for land higher than the Global average 
( 0.1% ) is to be expected and there is no inconsistency with 'our' value. 

Similar conclusions that most of the LCC variability comes from the subtropical 
oceans and is most likely due to TSI variations, causing changes in lower tropospheric 
static stability, have been made by Kristjansson et al. (2004). 

Another aspect of the Physics behind the correlation may be related to the rela- 
tionship between cloud height and cloud cover ( Cotton and Anthes, 1959 ). These 
workers estimate that changes of -0.3% in CC from a height of 5km ( mid MCC ) 
correspond to a change of +0.1°C at ground level ( i.e. cloud absorption of incoming 
radiation dominates ). In our case a change of +0.1°C at ground level corresponds 
to a change in LCC + MCC of -0.3%, i.e. the same result although using the whole 
of MCC is not really appropriate. Nevertheless, there is a similarity in the values. 
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5 Conclusion 



We advocate a scenario for the origin of correlations between CR and LCC, based 
on the parallel influence of solar activity. The solar irradiance rises with the sunspot 
number in the middle of the solar cycle. The radiation is strongest in the tropics 
and subtropics. Though the relative rise of the irradiance is small, and only about 
0.1%, it causes a rise of the mean surface temperature and an increase of the vertical 
convection flows of the heated air. The subsequent change in supersaturation of the 
air at different heights can cause the changes in LCC and MCC. Warm air from 
below 3 km rising to greater heights will cause the LCC to fall and MCC to rise. By 
this way the rise of convection flows leads to a considerable magnification (to ~2%) 
of the effect of enhanced solar irradiance. Formulating briefly, one can say that in the 
maxima of the solar cycles the updraft becomes stronger and this effect is strongest 
in the tropics and subtropics, as well as in the southern latitude bands where there 
is the largest fraction of area covered by the oceans. It is well known that the 
variations of solar activity are followed by the variations of CR intensity at Earth; 
the reduction of CR intensity coincident with the reduction of LCC is therefore by 
no means evidence of the causal connection between these two phenomena - they 
correlate with each other due to their common origin - the change of solar irradiance 
at the Earth. 
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Captions to figures 



Figure 1. The temporal behavior of CR (a), Global LCC (b) and the x 2 /naif value 
as a function of the time lag between CR and LCC curves (c). Dotted lines in (a) 
and (b) are the mean values and dashed lines are the 5-degree polynomial fits of the 
CR intensity and LCC respectively 

Figure 2. Correlation plots for variations of LCC and CR during the 22nd solar 
cycle: (a) long-term variations, (b) short-term variations. Dashed lines in both 
panels are linear regression lines, the full line in the (a)-panel is the best fit curve of 
Y = a + bX c type with coefficients a, b, c indicated inside the panel. The slope b of 
the regression line and the correlation coefficient r are indicated inside both panels. 

Figure 3. Left set of panels: the temporal behavior of MCC and LCC ( two upper 
panels ), correlation plot for their long-term and short-term variations ( two lower 
panels ). Right set of panels: the same as the left one, but for HCC and MCC 
respectively. 

Figure 4. Seasonal variation of the Global MCC (a) and LCC (b). 

Figure 5. The seasonal variation of the sea, land and Global temperature from 
( Global Surface Temperature Anomalies ). Numbers above the curves show the 
fraction of the total area occupied by the sea (71%) and land (29%). 

Figure 6. The temporal behavior of CR (a), LCC (b) and the correlation of their 
long-term variations during 22d and 23d cycles of the solar activity (c). 

Figure 7. The temporal behavior of the Global surface temperature (a), HCC (b), 
MCC (c) and LCC (d) during 22nd and 23rd cycles of solar activity. Dotted lines 
in (b), (c) and (d) panels show mean values of CC for 1985-2005 period of time. 
The dotted line in (a) panel shows the mean temperature during the last century 
1900-2000, which is equal to 13.86°C. It illustrates the higher temperature in the 
last two decades of this period - the so called Global Warming. Dashed lines show 
the linear fits of the temporal behavior of the surface temperature and CC. Full lines 
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are their 5-degree polynomial approximations. 

Figure 8. The latitude dependence of CC characteristics: (a) absolute values of LCC 
( open circles ), MCC ( full circles ) and HCC ( open stars ). (b) LCC, MCC and 
HCC correlations with CR ( Climax ). Notations are the same as in (a), (c) The 
sensitivity of MCC to LCC ( open circles ) and their correlation coefficient ( full 
circles ). 

Figure 9. Temporal behavior of the surface temperature in 1880-2005 years ( upper 
panel ) and its Fourier frequency spectrum of the global temperature for this time 
range ( lower panel ). The peak at 0.007 month -1 corresponds to ~1 1-year period 
coincident with 11.87 year solar cycle ( Sturrock, 2008 ). 

Figure 10. Temporal behavior of the LCC top pressure, fitted by linear ( dotted 
line ) and 5-degree polynomial ( full line ) approximations. 
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TEMPORAL VARIATION O E COSMIC RAYS AMD THE GLOBAL LOO IN THE CYCLE 22 




19 8 6 1988 1 990 1992 19 9 4- 199S 

Yea r 




— 12 — 10 — 8 — 6 — 4 — 2 O 2 4- 6 8 10 12 

Time difference between LCC and CR, months 



Fig. 1. 



21 



LCC vs. CR long — terr-n c o r r e I a t i o r 



b = O.T — 0.023 

r=0.53B 

best fit = O. 97827 




b = — 0.055H-/' — O . OS V 
r= — O . OQO 


(b) 







-O. 1 o 

OF? relative deviatior 



Fig. 2. 



22 
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SEASONAL VARIATION OF THE GLOBAL SURFACE TEMPERATURE 
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TEMPORAL VARIATION OF COSMIC RAYS AND THE GLOBAL LCC IN THE CYCLES 22 and 23 



O 
O 



1 1 1 1 1 1 1 1 1 1 

Neutro n counting rate at Climax ( a ) 


i i i 





1986 1988 1 990 1992 1994- 1996 1998 2000 2002 2004- 
Year 



' ' ' G'lobal Low Cloud Covir 




i i i i 


i i i 





1988 1988 1 990 1992 1994- 1998 1998 2000 2002 2004 

Year 

LCC vs. CR long — term correlation 



b=0.222 + /-0.0'34 
r=0.390 

— <& 


Co) 


'-v.. J "• 

° --'S'V ° ° °° 

. J = =r<? « 





J-O.2 



— O. 1 



CR relative deviation 



O. 1 



0.2 



Fig. 6. 



26 



TEMPORAL VARIATION OF THE GLOBAL TEMPERATURE, HOC, MCC AND LCC 
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LATITUDE DEPENDENCE OF CC CHARACTERISTICS 
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